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ABSTRACT 

The evolution of hot interstellar gas in cluster-centered cD galaxies and the inflow of gas from the 
surrounding galaxy clusters are strongly coupled. Cooling flows arise inside the cD galaxy because the 
deep stellar potential and stellar mass loss increase the gas density and decrease the radiative cooling 
time within the galaxy. Recent X-ray observations of M87 in the Virgo cluster and NGC 4874 in Coma 
reveal that the gas temperature beyond about 50 kpc from these cD galaxies is comparable to the virial 
temperature of the cluster, 3 or 9 keV respectively, but within the optical galaxy the temperature drops 
to the galactic virial temperature ~ 1 keV. We show that these steep thermal gradients on galactic scales 
follow naturally from the usual cooling inflow assumptions without recourse to thermal conductivity. 
However, most of the gas must radiatively cool ("dropout") before it flows to the galactic core, i.e. the 
gas must be multiphase. The temperature and density profiles observed in M87 and NGC 4874 can 
be matched with approximate gas dynamical models calculated over several Gyrs with either globally 
uniform or centrally concentrated multiphase mass dropout. Recent VMM observations of M87 indicate 
single phase flow at every radius with no apparent radiative cooling to low temperatures. Gas dynamical 
models can be made consistent with single phase flow for r > 10 kpc, but to avoid huge central masses 
of cooled gas, we assume that some distributed cooling dropout occurs near the center of the flow where 
the gas temperature is T ~ 1 keV. The evidence in X-ray spectra for multiphase cooling beginning at 
lower temperatures ~ 1 keV may be less apparent than for higher temperatures and may have escaped 
detection. However, even if the mass of cooled gas is distributed within r < 10 kpc, it is necessary 
that the mass of cooled gas not conflict with dynamical mass to light determinations. Because of small 
deviations from true steady state flow, we find that the standard decomposition methods used by X-ray 
observers to determine the mass flow M{r) may fail rather badly, particularly when the mass dropout 
decreases with radius in the flow. For this case the decomposition procedure gives the usual cooling flow 
result, M oc r, which is quite unlike the true variation of M{r) in our computed models. 

Subject headings: galaxies: elliptical and lenticular - galaxies: active - galaxies: cooling flows - 
galaxies: radio ~ X-rays: galaxies - X-rays: galaxy clusters 
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1. INTRODUCTION 

Both galaxy clusters and individual early-type galaxies 
contain virialized hot gas in which the radiative cooling 
time is much less than the Hubble time. According to the 
standard interpretation, as hot gas loses energy by emit- 
ting the X-radiation we observe, it slowly migrates inward. 
Since these "cooling flows" are nearly in hydrostatic equi- 
librium in the cluster or galactic potential, the tempera- 
ture and density profiles measured in the hot gas can be 
used to determine the underlying mass distribution. For 
this reason it is particularly interesting to study cooling 
flows in galaxies since we already have information about 
the mass distribution from optical studies within the half- 
light radius where the contribution of dark matter to 
the total mass is small. 

Our objective in the calculations discussed here is to de- 
termine the nature of the interaction between cluster and 
interstellar gas in bright cD galaxies centrally located in 
galaxy clusters. In our previous studies of X-ray luminous 
massive elliptical galaxies located in small galaxy groups 
(e.g. Brighcnti & Mathews 1999), we have assumed that 
the interstellar and circumgalactic gas have been relatively 
unperturbed since the galaxy first formed. By following 
the growth of cosmological perturbations, including star 
formation and feedback, we can compute the gas flow over 
most of the Hubble time, successfully reproducing the ra- 
dial dependence of hot gas density, temperature and abun- 
dance currently observed in massive group-centered ellip- 
tical galaxies. However, the hot gas surrounding cluster- 
centered elliptical galaxies has a much higher cluster virial 
temperature. As cluster gas flows into the central cD or 
elliptical galaxy, the interstellar gas is perturbed from that 
of more isolated or group-centered ellipticals, providing a 
usehil test of the many parameters involved in these calcu- 
lations. In particular, cluster-centered cD galaxies provide 
an excellent venue to explore the nature or reality of the 
mass "dropout" hypothesis, i.e. that the hot gas is multi- 
phase, having locally varying densities and cooling times, 
so that it cools ("drops out") throughout a large volume 
of the cooling flow. Unfortunately, the long term evolution 
of ambient cluster gas around cluster-centered cD galaxies 
is uncertain because of the likelihood of disruptive merg- 
ers in the cluster environment. Nevertheless, much can 
be learned if we assume that the cluster gas as currently 
observed near the cD galaxy has been relatively stable for 
times longer than several sound-crossing times, typically 
several Gyrs. 

In the traditional interpretation of cooling flows, it is 
assumed that the X-ray luminosity reflects energy lost 
from the hot gas which is not balanced by a heating source. 
As the radiating gas loses energy and slowly migrates to- 
ward the potential minimum, it is maintained near the 
local virial temperature by Pdv compression. In the sim- 
plest (single phase) cooling flows all the gas cools rapidly 
at very center of the E or cD galaxy where the gas density 
and radiative emission are greatest. This model clearly 
fails since strong localized central peaks of thermal X-ray 
emission are not generally observed and the total mass of 
cooled gas concentrated in the central regions would in- 
crease the stellar velocity dispersions considerably above 
those observed (e.g. Thomas 1986; White & Sarazin 1987; 
Mathews & Brighenti 2000). To alleviate these difficulties. 



it has been fashionable for many years to assume that the 
hot gas is "multiphase" and therefore cools throughout a 
larger volume, depositing its mass in a less concentrated 
manner (Nulsen 1986; Thomas et al. 1987; Allen et al. 
2001). 

Multiphase cooling can be incorporated into the equa- 
tion of continuity for cooling flows with a mass sink term: 

dp I d 2 s P /IN 

a^ + ^a^(-H = «*P*-?^- (1) 

The source term proportional to the stellar density rep- 
resents stellar mass ejected from evolving stars. The fol- 
lowing dropout or sink term simply removes gas from the 
flow at a rate characterized by the local radiative cooling 
time at constant pressure tcooi = 5mpkT/2^pA. This term 
is preceded by a dimensionless coefficient q{r) < 1 which is 
adjusted to best fit the observed X-ray surface brightness 
and other observations. We recognize the ad hoc character 
of this sink term, particularly in view of the well known 
faihirc of linear analysis to find instability in cooling flows 
(e.g. Balbus & Soker 1989). Nevertheless, the compelling 
astrophysical arguments against concentrated cooling at 
the center of the flow suggest either some form of dis- 
tributed mass dropout as in Equation (1) or an equally 
ad hoc heating effect that supplies thermal energy without 
driving a strong global wind. Ironically, we have recently 
shown (Brighenti & Mathews 2002) that nonlinear density 
perturbations and subsequent thermal instabilities occur 
naturally in cooling flows that are decelerated by heating 
or additional non-thermal pressure near the center of the 
flow. 

To estimate q we use the recent compilation of six well- 
observed cluster cooling flows from Allen et al. (2001) 
for which the local radiative cooling times are linear func- 
tions of radius, tcooi = tcooi,b{r /n), for r < n. Since 
rfc ~ 40 — 120 kpc is much greater than the central galac- 
tic half-light radius re ~ 10 kpc, a*/9* <C p/tcooi beyond 
the central galaxy. Multiplying the equation above by 
Airr'^dr, assuming quasi-steady state flow and substitut- 
ing p = bmpkT/2ntcooiK 



dM = dipuAirr^ 
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We see that tcooi k r implies M cx r as verified by Figure 
6 of Allen et al. who determine consistent mass flow or 
dropout profiles M(r) from both spectroscopic and sur- 
face brightness observations. Integrating this equation to 
r = rb and rearranging we find 



2M(r) ^in)tl,oi, 

207TmpkT rf 



(3) 



which is constant in r < in this simple approximation. 
Using data from the six well-observed clusters in Allen 
et al., we find {q) « 0.3 — 0.6. But the implication that 
tcooi oc r requires M oc r is very sensitive to the steady 
state assumption as we show in Section 4 below. 

Recently, XMM and Chandra have revealed new and 
unexpected attributes of galactic and cluster cooling flows 
(McNamara et al. 2000; Fabian et al. 2001; Fabian 2001). 
In several cluster cooling flows there is spectroscopic evi- 
dence for only partial cooling, down to ~ 2 keV, but no 
further (e. g. Peterson et al. 2001; Tamura et al. 2001; 
Kaastra et al. 2001; Molendi & Pizzolato 2001; Xu et al. 
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2001). This has led some to question the reahty of the 
g-term in the equation of continuity and to postulate in- 
stead some as yet poorly understood form of heating (e.g. 
Tucker & Rosner 1983; Binney & Tabor 1995; Tucker & 
David 1997; Ciotti & Ostrikcr 2000; Sokcr ct al. 2001). 
Although the traditional mass dropout hypothesis is be- 
ing challenged, at the present time there is no accepted 
alternative model with heating that can be compared in 
detail with observations. 

Nevertheless, some recent observations continue to sup- 
port the reality of traditional cooling flows. Buote (2000, 
2001) finds evidence of widespread X-ray absorption from 
warm gas at T ~ 10® K in many cooling flows including 
M87. FUSE observations of cluster A2597 have detected 
the OVI doublet from gas at T ~ 3 x 10^ K corresponding 
to a lower than expected M, but OVI lines are undetected 
in A1795 (Oegerle et al. 2001). Similar FUSE observations 
of elliptical galaxies find OVI emission consistent with (but 
possibly somewhat less than) the expected mass cooling 
rate M and the usual multiphase cooling flow model (e.g. 
Bregman et al. 2001). Perhaps some OVI line emission 
may have escaped detection because of broadening by tur- 
bulent velocities. Finally, Canizares et al. (1982) discov- 
ered line emission in M87 corresponding to temperatures 
of a few 10® K which they interpret at radiatively cooling 
gas. Although the evidence is somewhat contradictory at 
present, wo adopt provisionally the traditional assumption 
that cooling flows do in fact cool and seek values of q{r) 
for cluster-centered cD galaxies that agree best with recent 
X-ray observations. 

In the following we describe the transition flow as clus- 
ter gas enters M87 from the Virgo cluster and NGC 4874 
from the Coma cluster, both of which have been observed 
by XMM and/or Chandra. Remarkably, there is little evi- 
dence that the powerful, ongoing AGN in M87, with its jet 
and radio lobes, are heating the cooling flow at small galac- 
tic radii (e.g. Fabian 2001; Soker et al. 2001). Indeed, 
we demonstrate below that the gas temperature, density 
and entropy proflles in M87/Virgo and NGC 4874/Coma 
are approximately consistent with a normal galactic-scale 
cooling flows, requiring little or no additional local heat- 
ing. The interaction of cluster and galactic gas has been 
discussed previously by Thomas (1986), Meiksin (1990) 
and Bertin & Toniazzo (1995). 

2. m87 in VIRGO 

2.1. MSI/Virgo Mass Distribution 

Using ROSAT data, Nulscn & Bohringcr (1995) deter- 
mined the azimuthally averaged mass distribution of the 
Virgo cluster within 50' (247 kpc at distance c? = 17 Mpc) 
of M87. Their best-fitting model has a total density struc- 
ture 

poc[l + (r/a)2]-i (4) 

and mass profile 

Myirgo{r) = Mo[r/a + arctan(r/a)] (5) 

where a = 47.6 kpc and Mo = 5.90 x lO^^ 7^,/^ ^he 
mass distribution is sensitive to the observed gas temper- 
ature variation T(r) which for r > 85 kpc is approximately 
isothermal at ~ 2.8 keV (3.2 x 10^ K) (Nulsen & Bohringer 
1995). The observed X-ray surface brightness in Virgo is 



more consistent with this Myirgo{r) than with an NFW 
mass profile, although the diflference is not great. 

We determine the stellar mass distribution in M87 
(NGC 4486) from the B-band surface brightness /xb de- 
rived from the /iy and B-V photometry of Goudfrooij et 
al. (1994) within r — 132". Beyond this radius we use the 
V-band photometry of Schombert (1986) corrected to Hb 
with B-V = 0.80, and brightened by 5V = 0.5 to agree 
with the Goudfrooij et al. data in r ^ 132". The B-band 
surface brightness /is is transformed to surface mass den- 
sity using M/Lb = 8.2 (van der Marel 1991). (Note that 
Faber et al. (1997) find a much larger value, M/Lb = 23.) 
The stellar mass density, found by solving Abel's equation, 
can be fit with 

P*(r) = (r/rb)-^[l + (r/re)i-']-i[l + (r/rrf)]-°-4 Mq pc'^^ 

(6) 

where rt — 4.379, Vc — 1.513 and = 15 kpc. The stel- 
lar mass distribution M^^g? differs from a de Vaucouleurs 
proflle because of the unusually large core (break) radius 
in M87 and an outer cD envelope of moderate brightness. 
The potential is dominated by stars within 40 kpc. The 
density proflle in Equation (6) is in good agreement with 
that of van der Marel (1994). Finally, we assume that M87 
is surrounded by a dark galactic halo with an NFW pro- 
file and total mass Ms7,nfw = 1.0 x 10^^ Mq. While it is 
unclear if it is necessary to append a separate dark halo to 
the central cD, we note that dynamical models of cluster- 
centered E galaxies have been unable to reproduce their 
massive cD envelopes (e.g. Dubinski 1998). Fortunately, 
the mass Msi.nfw has only a minor influence on our re- 
sults because of the enormous mass of the Virgo cluster at 
every radius beyond the stellar half-light radius. 

2.2. Computational Procedure 

Each gas dynamical calculation begins with a determi- 
nation of the gas density profile of isothermal intracluster 
gas in hydrostatic equilibrium in the Virgo potential de- 
fined by Myirgo{r). At cosmic time ti the central galaxy 
and its dark halo are inserted into the Virgo environment. 
After this time stellar mass loss contributes to the local 
gas at a rate a^{t)p^{r) gm s~^ cm"*^. Here a^{t) = 
ct*,n{t/tn)~^'^ s~^, appropriate for a nearly coeval burst of 
star formation at high redshift where a*,„ = 4.7 x 10~^° 
s~^, and t„ = 13 Gyrs is the present age of the universe. 
Gas within the optical galaxy is heated by Type la su- 
pernovae which occur at a rate SNu(f) = SNu(i„)(i„/f) 
where SNu(f„) = 0.06 supernovae per 10^° Lb per 100 
years. These stellar mass loss and supernova rates are 
consistent with X-ray observations of massive E galaxies 
such as NGC 4472 that arc not centered in large clusters 
(e.g. Brighenti & Mathews 1999). 

After time ti, as the hot intracluster gas flows into the 
optical galaxy, it is cooled by radiation losses and mass 
ejected from evolving stars. Since mass-losing stars or- 
bit rapidly through the hot gas, we expect instabilities to 
thermally mix stellar and ambient gases. This assump- 
tion is confirmed by negative abundance gradients of a- 
clcmcnts observed in M87 (Molendi & Gastaldello 2001) 
which cannot have been produced by Type la supernovae. 
In large galaxies such as M87, the global potential is so 
deep that Type la heating has rather little effect on the 
thermal budget of the gas; the SNIa luminosity at 
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LsNia = 1-6 X lO'*'^ erg s~^, is considerably less than the 
total X-ray emission from within the optical image of 
M87. Following a brief transient (< 1 Gyr) just after M87 
is introduced into the Virgo environment, the gas flow set- 
tles into a subsonic, but secularly changing flow. The gas 
dynamical equations that we use arc identical to those dis- 
cussed in Brighenti & Mathews (1999) to which we refer 
the reader for further details. The equations are solved 
with a considerably modified ID spherical version of the 
ZEUS code (Stone & Norman 1992) using 360 logarithmi- 
cally increasing zones beginning with a central zone of size 
15 pc and extending to 2.1 Mpc. The appropriate initial 
temperature T{ti) for Virgo gas at time ti depends some- 
what on the nature of the assumed mass dropout profile 
q{r). If distributed mass dropout is small at large radii, i.e. 
g « 0, the isothermal temperature T{ti) should be close 
to that observed by Nulsen & Bohringer. But if w 1 at 
large radii, the temperature of the background flow T{ti) 
must exceed that observed, which is a weighted average 
over both background and cooling regions. 

The factors that influence the radial variation of gas 
density and temperature within the optical galaxy at time 
tn = 13 Gyrs arc: (1) the assumed mass dropout pro- 
file q{r) and therefore T{ti), (2) the central gas density 
n(r = 0, ti) of the hydrostatic cluster gas distribution just 
before the galaxy is inserted, (3) the time elapsed since this 
initial configuration t„— tj, (4) the specific stellar mass-loss 
rate and (5) thermal conductivity. As we show below, 
our computed values of n(r, tn) and T(r, tn) inside and 
near M87 are only weakly sensitive to the time tn — ti over 
which the gas fiow evolves after the initial setup, provided 
tn — ti is larger than ^ 3 Gyr, i.e. several sound crossing 
times. The time f„ — ti can be regarded as the (unknown) 
interval since the last major merger event in Virgo when 
the gas density and temperature may have changed ap- 
preciably. While it is unlikely that varies strongly with 
galactic radius, the normalization coefficient a*,„ can in- 
fluence n(r,tn) and T{r,tn) somewhat. It is unlikely that 
thermal conduction plays an important role in the thermal 
structure of M87/Virgo (e.g. Saito fc Shigeyama 1999). 
We find that thermal conduction at > 0.1 of the Spitzer 
(or saturated) rate cannot occur since the gas tempera- 
ture would be nearly isothermal at T « T{ti) ~ 3 keV 
even within ~ 30 kpc and this is not observed. Thermal 
conduction can help reduce the total mass of cooled gas in 
the central regions of M87, but this requires that the con- 
ductivity be reduced by a carefully chosen factor that pre- 
serves the observed temperature profile (e.g. Bertschinger 
& Meiksin 1986). Therefore, in our discussion below we 
concentrate mostly on the remaining factor that influences 
the; currently observed gas density and temperature distri- 
butions: the mass dropout coefficient q{r). 

2.3. Results 

In Figures 1 and 2 we compare the observed gas den- 
sity and temperature profiles in M87/ Virgo with results of 
our gas dynamical calculations. Unless otherwise noted, 
all calculations begin at cosmic time ti = 6 Gyrs and con- 
tinue to the present time t„ = 13 Gyrs. At time ti the 
hot cluster gas is in hydrostatic equilibrium in the mass 
distribution Myirgo{r). The cluster gas is isothermal at 
T{ti) = 3 keV (3.5 x 10^ K) for models with q{r) < 1 at 



r > 20 kpc and T{ti) = 3.4 keV (4.0 x 10^ K) for (globally 
multiphase) models with q = 1. The central gas density 
n{0,ti) = 0.003 cm~^ is chosen to agree with X-ray sur- 
face brightness observations at radii far beyond the optical 
extent of M87. 

If g = 1 throughout the flow, about half of the total 
X-radiation comes from the cooling elements and the rest 
comes from the background flow. We distinguish these two 
sources of emission in Figures 1 and 2 by plotting both 
the background density n^(r) (light solid line) and the 
density n{r) that would be determined from the combined 
emission from background and cooling regions (heavy solid 
line). This latter n(r) can be compared with the den- 
sity profile determined from observations. In the temper- 
ature plots we show three temperatures: the temperature 
Th{r) of the hot background gas flow at radius r (dotted 
lines), the mean local emission- weighted temperature T(r) 
including both background and cooling gas (light solid 
lines), and the emission- weighted temperature T{R) which 
is T(r) averaged along the line of sight, including emission 
from both background flow and cooling gas (heavy solid 
lines). R is the projected radius. If dT/dr > 0, T{R) is 
greater than T(r) for r = R because emission from hotter 
gas along the line of sight is included in the weighted av- 
erage T{R). The temperature T{R) reported by observers 
is not simply an emission-weighted mean along the line of 
sight, but is derived by fitting single temperature (IT) or 
composite tcmpcratm-e (2T or IT + CF) XSPEC mod- 
els to the projected X-ray spectrum. Nevertheless, in the 
following discussion we assume that our computed T{R) 
- which includes effects of dT/dr 7^ as well as cooling 
dropout when q - can be directly compared with ob- 
served temperatures. 

Over time f„ — f j = 13 — 6 = 7 Gyrs, cooling ffows do not 
form in Virgo without M87. If the evolution of the cluster 
gas after ti is calculated without M87, we find that T{R) at 
time tn rises from ^ 1.5 keV at the center, reaches 2.6 keV 
at r = 100 kpc, then slowly rises toward T{ti). The gas 
density is similar to n(r, ti) for r > 50 kpc, but within this 
radius it is much less than observed, n{r, tn) ^ 0.01 cm~^ 
and by time i„ no gas has cooled. However, when M87 is 
introduced at ti, the central gas density and X-ray surface 
brightness increase both because of the locally deepened 
(stellar) potential and because of the addition of stellar 
mass loss. We find that the deep potential of M87 is suffi- 
cient to create a cooling flow in M87/Virgo even if stellar 
mass loss is absent; this is possible since the Virgo virial 
temperature is only 2-3 times larger than that of M87. 
Stellar mass loss further cools the cluster gas and increases 
its density as it flows into the optical galaxy, lowering the 
radiative cooling time. 

2.3.1. Gas Flow with Constant and Variable q 

Figure 1 shows the resulting density and temperature 
profiles for the q = 1 model in which cooling dropout oc- 
curs uniformly throughout the flow. The computed den- 
sity variation n(r) (heavy solid line) agrees quite well with 
the gas density observed in M87/ Virgo. The density con- 
tinues to rise in r < 1 kpc, where the observational data 
are sparse, but in this region it is likely that the magnetic 
fleld is high enough ~ 100/xG to flatten the gas density 
profile. Owen, Eilek & Keel (1990) estimate a field of 
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~ 40/xG at r = 2 kpc. The projected model temperature 
profile T{R) also agrees reasonably well with observed val- 
ues but is somewhat lower for r 2 kpc and r > 45 kpc. 
M87 observations are more difficult near the center be- 
cause of non-thermal X-radiation from the active nucleus 
(e.g. Bohringer at al. 2001; Molendi & Gastaldello 2001), 
which may explain some of the temperature discrepancy 
at small r, along with a variety of other possible expla- 
nations (see below). The failure of the q = 1 model to 
match observed temperatures beyond ^ 45 kpc suggests 
that q may be less than 1 in this region since the observed 
temperatures in Figure 1 increase toward that of the un- 
cooled background gas T^(r) (dotted line) which is ~ 0.8 
keV hotter than T{R). 

To explore variable mass dropout, we show in Figure 2 
the results of a calculation with a cooling dropout that 
decreases with gas density and therefore galactic radius, 
q = exp[— 0.05/n(r)], but which is otherwise identical to 
the solution shown in Figure 1. As a function of radius, 
g « 1, 0.5 and 0.1 at r = 0.0, 2.8 and 7.6 kpc respec- 
tively. The agreement with observed density is very good 
and T{R) agrees adequately with observed values. It is in- 
teresting that the observed density n{r) and temperature 
T(i?) profiles can be fit for either constant or spatially 
variable q, particularly for r > 10 kpc. But the observed 
spectra of the two flows are quite different, as we discuss 
below. 

Inside ^ 2 kpc in Figure 1 or ^ 10 kpc in Figure 2 the 
computed temperature T{R) is somewhat lower than ob- 
served. It is unlikely that this discrepancy is due to shock 
heating by the M87 AGN or radio source. In our experi- 
ence the principal long term consequence of local heating 
is to lower the cc^iitral gas density with only a very small 
change in temperature. In any case, Chandra observa- 
tions provide little or no direct evidence of shock heating 
in radio galaxies (e.g. Fabian 2000). Several possible ex- 
planations for this temperature discrepancy are explored 
in the additional q = exp[— 0.05/n(r)] solutions shown in 
Figure 3. If the time since the initial M87/Virgo config- 
uration tn — U is ^ 4: Gyrs, we find that both T{R) and 
dT/dR are larger in R < 10 kpc. This is illustrated in 
Figure 3 with a solution att„ — = 13 — 10 = 3 Gyrs 
(long dashed lines) which is a poor fit to observed n(r) 
and T(r). However, for tn — ti > 4: Gyrs all T{R) profiles 
resemble that in Figure 2. Reducing the normalization of 
the stellar mass loss rate by a factor of 2 also has a 
beneficial, but limited affect on T{R) (short dashed line in 
Fig. 3). However, if the stellar mass to light ratio is twice 
that of van der Marel, M/Lb = 16.4 - and therefore in 
better agreement with the M/Lb of Faber et al. (1997) 
- the temperature profile is considerably improved (solid 
line). (Note: a* oc {M/Lb)~^ has been reduced by 2 in 
this model.) 

In Figure 2 T(R) is in excellent agreement with obser- 
vations in i? > 10 kpc. In this region the gas is essen- 
tially single phase, Tfi{r) — T{r), and the full line of sight 
variation ~ 2[T{R) - T{r)] « 0.4 - 0.7 keV is compara- 
ble to the cooling flow temperature range T^ax — Tmin in 
the XSPEC (Model C) analysis of Molendi & Pizzolato 
(2001) at r > 10 kpc. The multiphase appearance is due 
entirely to contributions of various single phase tempera- 
tures along the line of sight. At smaller radii, r 10 kpc, 



q increases, and T{R) is reduced by contributions from gas 
that is cooling below 1 keV. 

We have also repeated M87/Virgo calculations with no 
dropout at all, q = Q, so that all the gas cools right at 
the origin. The results of these calculations agree surpris- 
ingly well with observations at r > 1 kpc, particularly if 
a■^, is also lower. But there is a major problem. In r 1 
kpc both the gas density and temperature increase in the 
(point source) potential of the centrally cooled mass, pro- 
ducing a dense, luminous X-ray core. However, central 
masses ^10^° Mq and central peaks in T and are not 
typically observed, so cooling inflows without spatially dis- 
tributed dropout are unrealistic. 

2.3.2. Multiphase flow in M87/Virgo 

The detailed shape of the density and temperature pro- 
files in M87/ Virgo can be adequately reproduced with ei- 
ther constant or variable q, but this degeneracy can be 
broken by examining the X-ray spectra of these models in 
more detail. If g = 1 as in Figure 1, the line and continuum 
X-ray spectrum in r > 20 kpc should consist of two compo- 
nents of approximately equal flux: (1) a blend of emission 
from line of sight plasma at temperatures T/j ~ 2.8 — 3.4 
keV, and (2) emission from regions cooling at every radius 
through all temperatures T < T/j. However, if w as in 
our variable q model beyond 10 kpc, T{R) is determined 
by only the first of these two components. In the upper 
panel of Figure 4 we show approximate local spectra at 
r = 30 kpc for both constant and variable q models. For 
the q = exp[— 0.05/n(r)] fiow we use a single temperature 
(IT), solar abundance mekal model (using XSPEC version 
11.0.1) at T = 2.5 keV. For the q = 1 model spectrum we 
combine two bolometrically equal components: a 3 keV 
IT spectrum and a mekal cooling fiow (CF) spectrum with 
Tmax = 3 keV. The difference between the q = 1 and vari- 
able q spectra is most apparent near the Fe-L line complex 
near 1 keV which is broadened toward lower energies in the 
(IT + CF) spectrum for g = 1 flow. 

The current observational evidence for truly multiphase 
gas in M87/Virgo is contradictory. Buote (2001) reana- 
lyzed i?OS'^T PSPC data from M87/Virgo, including the 
important spectral range down to 0.2 keV, and concluded 
that a single temperature gas cannot explain the X-ray 
spectrum. Instead, Buote finds evidence of emission in 
the 0.2 - 0.4 keV band from warm (T 10^ K) gas in 
M87/Virgo at almost every radius ^ 70 kpc except, curi- 
ously, in the r = 5 — 10 kpc region where the multiphase 
flow appears to be suppressed. Buote also finds (oxygen 
edge) absorption in the 0.5 - 0.8 keV range throughout 
M87/Virgo consistent with absorption by the same warm 
gas that emits at lower energies. The mass of warm gas 
necessary for this absorption in M87 is enormous, com- 
parable to the total mass of hot gas observed, and it is 
unclear why the warm gas would not cool further (see also 
Bonamente et al. 2001). An important part of Buote's 
argument, shown in his Figure 4, is that if data below 
0.5 - 0.7 keV is neglected, as many authors have done, 
the absorbing foreground column Nh greatly exceeds the 
Galactic value; conversely, if the entire spectrum down to 
0.2 keV is considered, the foreground Nh is approximately 
Galactic at every projected radius in M87/ Virgo, but some 
warm gas absorption is required. Buote's ROSAT data 
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cannot easily discriminate between spectra of two temper- 
ature plasma and and a single temperature with cooling 
elements. 

This evidence for multiphase gas in M87/ Virgo is not 
supported by XMM data analyzed by Bohringcr ct al. 
(2001) which are spectroscopically superior to ROSAT. 
When they attempted a multiphase solution by combin- 
ing spectra from Th with that of a cooling gas, Bohringer 
ct al. found mass deposition rates of Af ^ 10 ±5 Mq yr"-"^ 
within 50 kpc and absorbing columns Nh ^ 5x 10^^ cm~^ 
that are about 30 times the Galactic value. 

Molendi & Pizzolato (2001) have reanalyzed 0.5 - 4.0 
keV XMM EPIC data from M87/Virgo and lend support 
to the conclusions of Bohringer et al.: (1) the only multi- 
phase attributes apparent in M87/ Virgo spectra are those 
introduced by the variation of single phase temperatures 
along the line of sight, and (2) there is no evidence for in- 
trinsic multiphase gas that is cooling locally. In their best 
multiphase cooling flow (Model B), Molendi & Pizzolato 
find (1) that the (spectroscopic) total mass deposition rate 
in M87/Virgo is less than 0.5 Mq yr~^, at least an order of 
magnitude less than that expected from L^, and (2) very 
large absorbing columns Nh- They suggest that earlier 
claims of large soft X-ray absorption are an artifact re- 
sulting from the multiphase assumption. But there remain 
problems with the model favored by Molcmdi & Pizzolato, 
Model C, in which the temperature range viewed along the 
sight is sufficient to account for the full temperature range 
observed at all R. They model this temperature range with 
a partial isobaric cooling flow from T^ax to T„j„ rather 
than by integrating over the known T(r) along the line of 
sight. But gas cannot just cool to some minimum temper- 
ature - ether locally as cooling dropout or during its radial 
flow - and accumulate there since that would eventually 
distort the entire spectrum toward this low temperature. 
This is the conundrum discussed by Fabian et al. (2001) 
who claim that "cooling flow" gas appears to cool to some 
intermediate temperature Tmin, then disappear. 

However, we must not overlook the Einstein FPCS 
observation of Fe XVII - Fe XXIV lines in the central 
- 2'(- 10 kpc) of M87 by Canizares et al. (1982) that 
reveal gas temperatures ranging from 3 x 10^ K down to 
2 — 4 X 10^ K. The volume emission measure from cooling 
gas suggests a cooling rate of 3 - 4 Mq yr~^. This is clear 
evidence in support of the multiphase model, at least in 
this central region. 

2.3.3. Spectroscopic Evidence for Variable q 

Returning to Figure 4, the easily visible low-energy ex- 
tension of the Fe-L line feature, which indicates cooling to 
very low temperatures, was evidently not seen at r ~ 30 
kpc by Molendi & Pizzolato (2001) in their XMM spec- 
tra of M87. This is in agreement with our variable q{r) 
model in which w in r > 10 kpc. However, in their 
(preferred) Model C analysis of XMM data in which the 
temperature cools through a flnite range, simulating the 
variation of T along the line of sight, Molendi & Pizzolato 
find temperature ranges of 1.0 < T < 1.8 keV at r = 3.9 
kpc and 1.4 < T < 1.9 at r = 7.5 kpc. These narrow 
ranges are inconsistent with intrinsic multiphase cooling 
in this region, with the observations of Canizares et al. 
1982, and with our variable q model. However, it may be 



more difficult to detect multiphase cooling when the gas 
cools from low temperatures, T ^ 1 keV, as in our variable 
q model within 10 kpc. 

To illustrate this, in the lower panel of Figure 4 we use 
XSPEC (ver. 11.0.1) to compare intrinsic 0.2 - 3 keV spec- 
tra at ^ 1 kpc, assuming solar abundances. The solid 
line shows the combined spectrum of a single tempera- 
ture (IT) gas at Th = 0.95 keV plus a cooling flow (CF) 
spectrum for T = 0.95 0.1 keV, with equal normal- 
izations. For comparison we show with a dotted line in 
Figure 4 the mekal spectrum of a single phase gas at the 
mean temperature of the IT -|- CF flow, T » 0.8 keV 
(cf. Fig. 2). The difference between these two spectra is 
much less than that between IT and (IT/i -I- CF) spec- 
tra at ^ 3 keV (upper panel) corresponding to i? = 30 
kpc where the Fe-L feature is widened by the cooling gas. 
To demonstrate conclusively that gas at ~ 1 keV is not 
cooling to very low temperatures, it would be necessary 
to be confident of the relative strengths of the Fe-L lines 
in 0.7 - 1.0 keV which may not be possible at present. 
Alternatively, the line ratio OVII(0.574 keV)/OVIII(0.652 
keV) increases with decreasing temperature below 1 keV 
(Landini & Monsignori Fossi 1990), which is apparent in 
Figure 4. However, OVII(0.574 keV)/OVIII(0.652 keV) 
could be artificially increased by oxygen-edge absorption 
(Boute 2001), which is a further possible complication. 
On balance, therefore, these spectroscopic considerations 
favor the variable q{r) model for M87/Virgo, in which q 
becomes effectively nonzero only within 10 kpc where a 
multiphase spectrum would be expected, but which may 
be difficult to observe. 

2.3.4. Where is the Cooled Gas? 

A potentially serious problem with the variable q model 
is that the centrally concentrated mass dropout may con- 
flict with the stellar mass to light ratio of M87. In his 
global dynamic models for M87, van der Marel (1991) 
found M/Lb = 8.2 while Faber et al. (1997) find 
M/Lb = 23 using the King method applied to the galactic 
core (Richstone & Tremaine 1986; Tremaine et al. 1994). 
In these estimates M/Lb is always assumed to be uni- 
form with radius - which would not be expected if mass 
dropout is important. The discrepancy in M/Lb could 
indicate that dark dropout mass exists mostly within the 
M87 core, although van der Marel (1994) also finds a stel- 
lar M/Lb = 8.2 in stars near the central black hole and 
Merritt & Oh (1997) find M/Lb = H- The black hole 
mass considered by van der Marel, M^h = 3 x 10^, is only 
three times less than the total stellar mass Mf, = 1.0 x 10^*^ 
Mq within the core (or break) radius = 626 pc (7.61"). 
This black hole mass is confirmed by observations of a thin 
HII disk in Keplerian rotation to within 4 pc of the center 
(Marconi et al. 1997; Macchetto et al. 1997). 

Most of the dark matter deposited in the M87 core by 
the variable q solution occurs at very early times. To study 
this in more detail we use an "alternative" M87/ Virgo 
model in which the initial evolution of M87 occurs outside 
Virgo and the galactic scale cooling flow deposits mass ac- 
cording to q — exp(— 0.05/n). By time ti — 6 Gyrs when 
M87 is inserted in Virgo, the total mass of cooled (dropped 
out) gas is already M^o = 2.3 x 10^° Mq or 2.3Mb. At 
these early times our schematic model for galaxy evolution 
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is very uncertain. Since the stars in M87 have significant 
metal enrichment, some of the enriched galactic hot gas 
that cooled at early times must have formed into second 
generations of old luminous stars observed today (Nulsen 
& Fabian 1997). But it is miclear how long cooled gas 
formed into optically luminous stars. Many authors have 
suggested that at the present time only low mass, optically 
dark stars are forming in cooling flows, but other lines of 
evidence, such as the stellar H/? index, suggest that lumi- 
nous stars continue to form at the present time (Mathews 
& Brighenti 1999). Furthermore, the precise value of the 
total mass of cooled gas at time U depends critically on the 
unknown cooling and merger history of M87 before ti and 
whether or not q = exp(— 0.05/n) applies at these early 
times, as we have assumed. By time i„ = 13 Gyrs, in the 
variable q flow 5.3 x 10^° Mq has cooled. 

If the cooling dropout hypothesis is correct, the cur- 
rently observed stellar mass to light ratio must reflect in 
part a contribution from cooled gas. If luminous stars 
form continuously from cooled gas, the mass to light ra- 
tio could in principle be similar to that of the old stellar 
population, although an exact match seems unlikely. On 
the other hand, if only non-luminous stars formed in our 
q = exp(— 0.05/n) flow, then at time ti = 6 Gyrs M/Lb 
would have already increased to 18 within the break radius 
rb, ignoring passive cosmological luminosity evolution. By 
time tn = 13 Gyrs, M/Lb within increases further to 28 
which is unacceptably high. Notice, however, that in our 
models the mass dropout continues to increase with ris- 
ing gas density in r ^ where about half of the dropout 
occurs for the variable q flow. It is more likely that the 
gas density in M87 levels off in r < 1 kpc similar to other 
large elliptical galaxies in Virgo (e.g. Brighenti & Math- 
ews 1997), possibly due to large central magnetic fields or 
the central AGN. If so, mass dropout within would be 
much reduced. More generally, the success of the variable 
q flow does not depend critically on the exact profile of q(r) 
or Mdoif) within ~ 7 kpc. To consider the most optimistic 
case, suppose all the cooled gas at time tn were distributed 
uniformly with old stars, Mdo{r) oc M*(r), within 7 kpc, 
then the central M/Lb would increase from 8.2 to only 
8.6. In order to accommodate any increase in M/Lb, it 
is necessary to assume that the underlying old stars in 
M87 have M/Lb less than our adopted 8.2 and that the 
dark dropout has increased it to this value. Finally, mass 
dropout in the q ~ I flow is more spatially diffuse: for 
optically dark dropout the mass to light ratio within rj, is 
M/Lb = 12 at ti and M/Lb = 14 at i„ = 13 Gyrs. 

It is well known that thermal conductivity can reduce 
the cooling rate in the centers of central cD galaxies in 
galaxy clusters, but the thermal conductivity must be 
fine-tuned to lower M without unrealistically flattening 
the temperature gradient (Bertschinger & Meiksin 1986; 
Meiksin 1988; Bregman & David 1988; Rosner & Tucker 
1989). We have repeated q = 1 and variable q{r) calcula- 
tions for M87/Virgo, assuming the "alternative" evolution 
discussed above and including thermal conduction. We as- 
sume that the Spitzer or saturated thermal conductivity 
is reduced uniformly by a factor / to simulate suppression 
by tangled magnetic flelds. In Figure 5 we show that the 
density and temperature profiles in M87/Virgo can be fit 
reasonably well with / = 0.02 for g = 1 flow and / = 0.08 



for variable q flow. If / is less than these values, the con- 
ductivity has little effect and the final flows resemble those 
described earlier without conduction. 

Conductive heating reduces cooling dropout. For the 
variable q solution with / = 0.08 the total mass of cooled 
gas is 3.1 X 10^" Mq, about half that without conductivity. 
Unfortunately, these results are very sensitive to the factor 
/. The dashed lines in Figure 5 show that if / is changed 
to 0.05 {q = 1) or 0.15 (variable q), the gas temperature is 
unacceptably high near the center. Since there is no obvi- 
ous reason why the conductivity should be lower by just 
the factor / required to reduce \M\, these models cannot 
be accepted. 

3. NGC 4874 IN COMA 

Evidence that the Coma cluster is not in dynamical 
equilibrium was noticed at ROSAT resolution by Briel, 
Henry & Bohringer (1992) and Briel & Henry (1997) who 
found significant subclustering and asymmetric temper- 
ature variations. Direct ROSAT images (Dow & White 
1995) and the wavelet analysis of X-ray data by Biviano 
et al. (1996) clearly show that the center of Coma con- 
sists of two subclusters, each centered on the bright cD 
elliptical galaxies NGC 4489 and NGC 4874, separated in 
projection by 7.1' . Recent XMM data confirm the exis- 
tence of enhanced, ^ 100 kpc X-ray emission around these 
two galaxies (Arnaud et al. 2001; Briel et al. 2001) and so 
do the Chandra images of Vikhlinin et al. (2001). We as- 
sume here that the two bright E galaxies are each centers 
of merging subclusters so the ambient dark matter and 
gas peaked around each galaxy is locally comoving with 
the galaxies as they orbit in Coma, i.e. there is no local 
ram pressure effect near the cD galaxies. The transient, 
dynamically disturbed nature of Coma is consistent with 
the observation that Coma is among the minority of clus- 
ters that do not have the characteristic centrally peaked 
X-ray cooling flow structure. In view of the unrelaxed na- 
ture of Coma, we begin our gas dynamical calculations at 
a more recent cosmic time ti = 10 Gyrs, but tn — ti is still 
longer than the relevant sound crossing times. 

The gas temperature within a few kpc of NGC 4489 and 
NGC 4874 decreases dramatically from the virial temper- 
ature of the Coma cluster, T « 8 — 9 keV at r > 7 kpc 
to the virial temperature of the embedded galaxy T ~ 1 
keV at r < 4 kpc. (Arnaud et al. 2001; Vikhlinin et al. 
2001). Because of its larger temperature excursion, Coma 
is a more extreme example of the interaction between clus- 
ter and galactic gas than M87/ Virgo. Yet, as we show 
below, the steep temperature gradient observed can also 
be understood with the classical cooling flow hypothesis. 
According to the Chandra observations of Vikhlinin et al 
(2001), the emission- weighted temperature near NGC 4874 
should increase from 

(Ti) = 1.0 ± 0.04 keV inside 7" (3.40 kpc) (7) 

to 

(T2) = 8.9±0.4keV in 10" -60" (4.86-29.2 kpc). (8) 

(T2) is a line of sight average. In deriving (Ti), Vikhlinin 
et al. removed background emission from 15 - 30", but 
not the local background in the 7 - 15"region. Neverthe- 
less, (Ti) approximately represents an average tempera- 
ture over the physical radius interior to 7" or 3.4 kpc. We 
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assume a distance to Coma of 106 Mpc based on Hq = 70 
km s"^ per Mpc and v = 7461 km (Faber et al. 1997). 
At this distance 1"= 0.486 kpc. 

3.1. NGC 4874/Coma Mass Distribution 

The two bright elhptical galaxies in Coma have very 
similar optical and X-ray properties. However, we con- 
centrate here on NGC 4874 because it is more ccintrally 
located (in projection) in the cluster and the surround- 
ing X-ray isophotes are more symmetrical. We assume 
that the mass distribution is also spherically symmetric 
about NGC 4874. Since our interest is in the gas flow 
pattern within ~ 100 kpc from the center of NGC 4874, 
it is not necessary to approximate the complete mass dis- 
tribution of Coma with high precision. We adopt a large 
scale isothermal total mass profile Mcoma{f) as given in 
Equation (5) with Mo — 2.8 x lO^'' Mq and core radius 
a = 300 kpc, which is atypically large due to dynamically 
transient substructure. This mass confines the cluster gas 
with an X-ray surface brightness distribution similar to 
that observed by Bricl et al. (1992) out to r ^ 1 Mpc 
and beyond. In addition, we adopt an NFW dark matter 
halo for the galaxy, M4874^„/^ = 1 x 10^^ Mq, although 
this mass has almost no effect on our solutions since the 
galactic virial temperature is so much less than that of the 
Coma cluster. 

The stellar mass distribution is found from the R-band 
surface brightness measurements of Peletier et al. (1990), 
which vary linearly with the log of the projected radius R, 

= 17 + 3.125 log i? mag/D". (9) 

The corresponding B-band brightness /xb is found assum- 
ing B — R = 1.67 typical for E galaxies (Fukugita, Shi- 
masaku & Ichikawa 1995). As before we convert /is into 
physical units and invert with Abel's equation to find the 
following approximate stellar space density 

p*(r)=p4r6)(r/r6)-0-^« r<n (10) 

= p^{rb){r/rby^-^'^ n < r < Vmax (H) 

where p^n) = 2.137 Mq pc'^, n = 2.67" (1.30 kpc) is 
the break radius (Faber et al. 1997) and Vmax = 33.71r(,. 
This />, (r) is very close to the stellar mass density shown in 
Figure 4 of Harris et al (2000) and we have used M/Lb = 8 
from their paper. The stellar mass distribution of NGC 
4874 dominates the gravitational potential in r ^ 12 kpc. 

3.2. Results 

As before, we begin our NGC4874/Coma calculation by 
establishing an isothermal gas background in the cluster 
with temperature Tcoma ~ 8 — 9 keV. Initially the gas is 
in hydrostatic equilibrium in the approximate dark halo 
of Coma, Mcomair)- At time ti = 10 Gyrs we place 
NGC 4874 and its dark halo at the center of the cluster, 
and begin the gas dynamical calculation with the same 
stellar mass loss and supernova parameters that we used 
for M87/Virgo. We adopt a more recent starting time 
ti = 10 Gyrs for NGC4874/Coma because Coma is dy- 
namically unrelaxed compared to M87/Virgo. The tem- 
perature Tcoma and central density n(0, ti) at time ti de- 
termine the pressure of the cluster gas near NGC 4874 
and are adjusted to agree with the large scale gas density 
distribution in Coma. 



For a fully multiphase flow model with = 1, we use 
n(0,t,) = 0.0015 cm-3 and T^oma = 13 keV (1.5 x 10^ K). 
Wh(ui (7 = 1, the; background cluster temperature must 
be significantly higher than temperatures observed since 
the latter also includes cooling gas at every radius. Figure 
6 shows the gas electron density and temperature for the 
(7=1 model after it has evolved to time t„ = 13 Gyrs. 
In discussing their Chandra observations of NGC 4874, 
Vikhlinin et al. (2001) show the X-ray surface brightness 
distribution E^;, but not the density profile n{r). Since 
NGC 4874 was 2.8' from the center of the Chandra field, 
the spatial resolution is degraded to a half-energy radius of 
~ 0.9". Figure 7 shows how the X-ray surface brightness 
Hx (0.2-5 keV) of our q = 1 solution (solid line) is altered 
by such a Gaussian point spread function (dashed line) us- 
ing the integration procedure of Pritchet & Kline (1981). 
In Figure 7 we plot the (unnormalized) surface brightness 
observations of Vikhlinin et al. (2001) (circles) and Dow 
& White (1995) (squares) by adjusting the data vertically 
for the best fit. The fit is acceptable except perhaps near 
i? ~ 2 — 3 kpc where the model emissivity is too high, but 
this is the region where Chandra observations become less 
accurate in 0.5 - 2 keV. 

The physical T{r) and projected T{R) temperature pro- 
files for the g = 1 calculation shown in Figure 6 can be used 
to determine the mean temperatures in the two regions in 
Equations (7) and (8): (Ti) = 1.0 keV and (Ta) = 7.8 keV. 
The latter temperature is less than the value {T2) ~ 8.9 
keV observed by Vikhlinin et al. (2001), but is in better 
agreement with temperatures from Arnaud et al. (2001) in 
this same region. (T2) is also much less than the assumed 
temperature in the background gas, Tcoma = 13 keV. Also 
shown as filled circles in Figure 6 are local tcnnperatures 
determined by Vikhlinin et al. which agree rather well 
with the (unprotected) T(r) in this model. 

We also consider an identical flow, but with spatially 
variable dropout, q = exp[— 0.05/n(r)], and initial density 
n(0,ti) = 0.003 cm and temperature T^oma = 9 keV. 
These solutions shown in Figure 8 become essentially sin- 
gle phase {q « 0) for r > 5 kpc. Beyond a few kpc the 
apparent temperature T{R) exceeds the temperature at 
physical radius r = R, T(r), because of line of sight projec- 
tion. The mean temperature in the two regions as in Equa- 
tions (7) and (8) are: (Ti) = 1.0 keV and (T2) = 8.5 keV. 
(72) is only 0.4 keV less than that observed by Vikhlinin et 
al., within the margin of observational error. The surface 
brightness 'Sx (0-2 - 5) keV of this variable q model shown 
in Figure 7 is also a better fit to the observations at all R. 

If the dark and stellar mass of NGC 4874 are ignored 
as well as stellar mass loss and supernovae, we find that 
no galactic-scale cooling flow develops in our Coma model 
even after computing for tn — ti = 13 Gyrs when the cen- 
tral temperature has dropped by only ~ 0.5 keV below 
its initial value. Although there is no cluster-scale cooling 
flow in Coma, more relaxed clusters with Tyi^ ~ 9 have 
higher central densities and can develop cooling flows even 
without a central cD. 

We have also performed NGC 4874/Coma calculations 
including thermal conductivity as suggested by Vikhlinin 
et al. (2001). Conduction with unmodified Spitzer and 
saturated coefficients can definitely be ruled out since the 
gas temperature remains isothermal at ~ 8 keV through- 
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out NGC 4874 even for large t„ — ti. Our results do not 
support the suggestion of Vikhlinin et al. that the sharp 
drop in gas temperature from 8 keV to 1 keV over a few kpc 
can be understood simply by setting the conductive heat 
flow (using a conductivity at the classical Spitzcr level) 
at some radius equal to the interior radiative X-ray lumi- 
nosity near NGC 4874. Instead we find that the coeffi- 
cients for saturated and unsaturated conduction must be 
reduced by more than 300 before (Ti) drops below 1.5 keV. 
As the conductivity is reduced further, our solutions ap- 
proach the observed (Ti) = 1 keV just as the conductivity 
term in the thermal energy equation becomes too small 
to influence the temperature profile at all. This result 
for conduction in NGC 4874/Coma is essentially the same 
for either constant or variable q flows. We conclude that 
thermal conduction in NGC 4874/Coma is magnetically 
suppressed to levels of negligible relevance. 

In summary, Chandra and XMM temperature profiles 
near NGC 4874 can be fit reasonably well with either con- 
stant or variable dropout q, although variable q flows are 
in better agreement with Chandra observations. As with 
M87/Virgo, the preference between these two models must 
be based on the evidence in high quality X-ray spectrum 
for or against multiphase flow. 

4. VARIATION OF m(r) 

As we mentioned in the Introduction, the observed den- 
sity and temperature profiles can be used to estimate 
tcooi{f) and therefore the integrated mass flow M{r) at 
radius r for steady state cooling flows. Typically, it is 
found that M oc r for r < ~ 40 — 120 kpc (e.g. Allen 
et al 2001). The decomposition procedure often used by 
X-ray astronomers to flnd M{r) from the observed T(r) 
and n(r) is based on dividing the observed cluster into 
an ensemble of concentric spherical shells and determining 
the bolometric X-ray emission from each shell, assuming 
steady state inflow with no stellar mass loss or heating 
(e.g. Arnaud 1988; Allen et al. 2001). The luminosity of 
the jth shell is 



L,=AM,(If,+A$,) 



A$,) 



(12) 



where AMj is the mass deposited in the jth shell, Hj = 
5kTj/2ij.mp is the enthalpy of the jth shell and A$j is the 
gravitational energy released in crossing the shell. The 
physical signiflcancc of each term is explained in Arnaud 
(1988). It is interesting to compare the mass flow 

Afi2(r,) = 5^AM, 

1=1 

found from Equation (12) with M(r) from our computed 
models. 

The heavy solid lines in the upper two panels of Figure 
9 show M(r) = pu^wr^ for our q = 1 and variable q flows 
in M87/Virgo at time t„ ~ 13 Gyrs. Also shown are the 
rates Mdo{r) = Jq ^''Tr'^{qp/tcooi)dr that mass is dropping 
out of the flow (dotted lines); note that this corresponds to 
the integrated spectroscopic mass deposition rate. If the 
flow were perfectly steady state M(r) and Mdo{r) should 
be equal as they nearly are for the q = 1 solution. For 



variable q flow, however, g « beyond ^5 — 10 kpc, so 
the integrated dropout mass levels off (dotted line in cen- 
ter panel). For this flow M(r) w M^gir) only holds for 
r ^ 20 kpc and beyond this radius they diverge in mag- 
nitude and sign because of small deviations from steady 
flow. The approach to steady flow is accelerated m q = 1 
flow (upper panel) because mass dropout tends to damp 
small transient velocities, i.e. dropout makes subsonic flow 
even more subsonic. 

In Figure 9 we also illustrate the mass deposition rate 
M-i2{t) determined by applying Equation (12) above to our 
computed flow at time tn (dashed lines). In solving this 
equation, we used the same gravitational potential used in 
our M87 flow calculations and evaluated the enthalpy us- 
ing the temperature Th{r) of the background flow at time 
tn- For (? = 1 flow Mi2(r) agrees quite well with both 
M(r) (solid line) and Mdo{r) (dotted line) in spite of the 
somewhat transient nature of the flow and the presence 
of stellar mass loss in r ^ 10 kpc. However, Equation 
(12) fails to capture either the mass flow or mass deposi- 
tion profiles in the variable q flow (Fig. 9 center panel). 
Even when Equation (12) is inaccurate, we seem to find 
the classic result, M12 oc r. In the variable q solution 
beyond about 10 kpc, small undamped radial motions re- 
main at time tn which have almost no infiuence on n{r) or 
T(r), but which radically alter M(r). We conclude that 
Equation (12) is not a reliable means of estimating M{r) 
in clusters that have small deviations from steady flow. 
Indeed, this may apply to all clusters observed. 

For comparison we also plot (lower panel of Figure 9) 
M(r) and Mdoir) based on a variable q calculation in 
which the interstellar flow in M87 is computed as an iso- 
lated galaxy (as in Brighenti & Mathews 1999), then in- 
serted into Virgo at time ti = 6 Gyrs, i.e. the "alternative" 
evolutionary model previously discussed. This computa- 
tional procedure introduces fewer velocity transients. The 
resulting model generates n(r) and T(r) proflles at time 
tn that are virtually identical to those plotted in Figure 2, 
but the corresponding M(r), shown as a solid line in the 
lower panel of Figure 9, is radically different from M{r) in 
the center panel. M{r) in the alternative solution (solid 
line) lies closer to Mi2(r) (dashed line) but both deviate 
from the true mass dropout profile Mdoix) (dotted line). 
Asymptotically for large tn — ti, we expect both M{r) and 
Mi2(r) to approach Mdo{r), but this may take longer than 
a Hubble time. The approximate agreement between M(r) 
and Mi2('') estimated from Equation (12) in r ^ 50 may 
be fortuitous. Equation (12) is nonetheless incapable of ac- 
curately detecting the presence of radially dependent q{r) 
in our models. 

Finally, in Figure 10 we show M, Mdo and M12 for the 
two models with thermal conductivity suppression factors 
/ chosen to best fit the n{r) and T(r) data in Figure 5. 
Both of these are alternative models in which the flow in 
M87 is computed before inserting into Virgo. When q = 1 
the conductivity has little effect on the various mass flow 
rates. But flows with q = exp(— 0.05/n) show that Mdo 
can be reduced by about 3 from the corresponding calcula- 
tion without thermal conduction (lower panel in Figure 9). 
As before, Equation (12) is inaccurate when applied to so- 
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lutions with variable q and continues to indicate M12 oc r, 
predicting mass deposition rates that are much too large. 

5. CONCLUSIONS 

Recent Chandra and XMM observations of M87 in the 
Virgo cluster and NGC 4874 in Coma show that the gas 
temperature decreases dramatically from the cluster virial 
temperature (3 or 9 keV respectively) at ~ 50 kpc down 
to the galactic virial temperature 1 kcV) on scales of 
a few kiloparsccs. We have performed a series of gas dy- 
namical calculations to determine the physical nature of 
this sharp thermal gradient. 

In addition, recent observations have presented some- 
what contradictory information about the important ques- 
tion of single phase or multiphase cooling flows. In multi- 
phase flows some of the hot gas cools throughout a large 
volume to avoid accumulation of large masses of cooled 
gas right at the center of the flow. Multiphase flow can 
be identified from spectral observations that show a large 
range of temperatures at a particular projected radius in 
the flow, larger than that expected from radial temper- 
ature variations projected along the line of sight. As a 
result of local cooling, gas mass "drops out" from multi- 
phase flows and the gas density in the flow decreases as 
{dp/dt)do = —Qp/tcooi where tcooi is the local radiative 
cooling time and < q{r) ;^ 1 is a mass dropout function 
that is chosen for a best fit to the observations. 

Our gas dynamical studies of the interaction of cluster 
and interstellar gas in cluster-centered cD galaxies are nec- 
essarily approximate because of unknown stochastic merg- 
ing events that may have influenced the cD environment 
in the past. Nevertheless, it is possible to estimate the 
flow of cluster gas into central cD galaxies provided the 
present day properties of the cluster gas have remained 
approximately constant over times for the flows to become 
quasi-steady on galactic scales. In most of the gas (cool- 
ing) flows calculated here we abruptly place the cD galaxy 
into the cluster at some time ti after which there is a tran- 
sient flow that relaxes to the desired solution after several 
Gyrs. We have also performed "alternative" evolutionary 
calculations in which the central cD galaxy evolves by it- 
self from early times, similar to our galaxy group model 
for NGC 4472 (Brighenti & Mathews 1999), and is placed 
in the cluster environment at time ti. The results of these 
two types of flow calculations are very similar except the 
total mass of cooled gas is somewhat larger in the alterna- 
tive flows. 

As a result of our calculations for M87/ Virgo and 
NGC4874/Coma, we reach the following conclusions: 

(1) The radial variation of gas temperature and density 

on scales of the cD galaxy can be understood from the 
usual cooling flow assumptions (e.g. stellar mass loss, spa- 
tially distributed mass dropout, etc.) without requiring 
additional heating sources such as thermal conduction. 

(2) If we perform our gas dynamical calculations for the 
Virgo cluster, assuming isothermal gas at T ~ 3 keV with 
the currently observed large scale density distribution, we 
find that a fully mature, large scale cooling flow does not 
develop without the presence of the central cD galaxy. Af- 
ter many Gyrs of calculation, a cD-less single phase flow 
in the Virgo potential exhibits line of sight averaged gas 



temperature profiles T{R) with a broad central minimum, 
but no gas cools. However, a fully complete cosmological 
calculation of the Virgo cluster is expected to produce a 
cooling flow with or without a central galaxy (e.g. Sugi- 
nohara & Ostriker 1998). 

(3) In our model for the richer Coma cluster with virial gas 
temperature ^ 9 keV, we find that the central density is 
sufficiently low that no cluster scale cooling fiow develops 
in a Hubble time without the central cD. In the presence 
of the cD galaxy, a cooling flow forms on a galactic scale, 
just a few kpc in size. This conclusion is probably not 
general. Since Coma is dynamically disturbed, its current 
central cluster gas density is atypically low. Relaxed clus- 
ters of similar richness have higher central gas densities 
and may form global cooling flows even without a central 
cD galaxy. 

(4) The radial gas density and temperature profiles in both 
M87/Virgo or NGC4874/Coma can be understood either 
in terms of models with large constant mass dropout at 
every radius {q = 1; multiphase flow), or with a variable 
dropout coefficient q{r) that becomes very small (single 
phase flow) beyond ^ 10 kpc. High quality X-ray spectra 
can detect the full range of local temperatures in the flow 
and thus determine the single or multiphase character of 
the flow. An observation of the oxygen ratio O VII (0.5 74 
keV)/OVIII(0.652 keV) would be useful to detect gas cool- 
ing below 10'' K. However, currently available XMM obser- 
vations of M87 have been unable to detect unambiguous 
evidence for multiphase cooling, even in the central regions 
of M87. 

(5) If we insist that the gas must cool somewhere in the 
fiow, we suggest that this occurs within the optical ex- 
tent of the central cD galaxy where the gas tempera- 
ture is close to the galactic virial temperature ~ 1 keV. 
At such low temperatures it is observationally more diffi- 
cult to detect intrinsic multiphase flow than in hotter gas 
more distant from the cD. These remarks are specific to 
NGC4874/Coma and M87/Virgo; in more massive cooling 
flow clusters, the cooled dropout mass may be distributed 
over a larger spatial region. 

(6) Unfortunately, when most of the radiative cooling is 
restricted to the central regions of the cD galaxy (r ^ 10 
kpc), the mass of cooled gas there may be inconsistent 
with mass to light ratios determined from stellar kinemat- 
ics. At early times when most of the cooling mass dropout 
occurred, the cooled gas probably formed into luminous 
stars. Stellar H/3 observations suggest that luminous stars 
may continue to form from cooling dropout. If the dropout 
mass is nonluminous, as many authors assume, then it is 
easy to violate observations of the stellar mass to light ra- 
tio unless the nonluminous mass is widely dispersed over 
the inner ~ 7 kpc in M87. 

(7) In this paper we have sought ways to hide the multi- 
phase nature of gas cooling within 10 kpc of the center of 
M87 and to hide the presence of large masses of optically 
dark cooled gas there. In spite or our limited success, it 
is difficult to fully embrace this approach. If M87 were in 
a much richer cooling fiow cluster, it may be impossible 
to assert that all of the cooled gas accumulates within the 
central 10 kpc. 
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(8) Wc confirm early studies that the coefficient for ther- 
mal condnctivity must be unrealistically fine-tuned or en- 
tirely negligible to lower the mass deposition rate without 
disagreeing with the observed temperature profiles. 

(9) For the cooling flow models considered here we find 

that the mass flow rate M{r) = puAur'^ can be very sen- 
sitive to small transient velocities in the flow that have 
no consequence for n{r) and T(r). Furthermore, the often 
used Equation (12) is incapable of accurately detecting the 
presence of radially dependent mass dropout q{r) in flows 
that are not completely relaxed, which may apply to the 
majority of cooling flows. In particular, when flows be- 
come nearly single phase at large radius, as indicated by 
recent XMAi data, our models suggest that Mi2(r) found 
from Equation (12) predicts a larger mass dropout at the 



current time than really exists in the flow. 
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Fig. 1. — Gas density and temperature for the q = 1 model of M87/Virgo at time t„. Upper panel: Density profiles rt(r). Light solid 
line: density of the background hot gas, heavy solid line: mean local density including emission from cooling regions. The data are taken 
from Nulsen & Bohringer (1995). Lower panel: Temperature profiles for the 9=1 model, dotted line: Thir), the background, uncooled gas; 
light solid line: T{r), the mean emission- weighted local temperature at physical radius r including cooling regions; heavy solid line: T{R) 
emission-weighted mean of T(r) along the line of sight at projected radius R = r. Data are from Molendi Sz Gastaldello (2001) {filled circles) 
and Bohringer et al. (2001) {open circles). 

Fig. 2. — Gas density and temperature for the variable q model of M87/Virgo at time tn- Line types and data are identical to those in 
Figure 1. 

Fig. 3. — Gas density and temperature in M87/Virgo at time tn for three variants of the variable q model. In the lower panel only the 
projected mean temperature T{R) is plotted. Data are identical to that in Fig. 1 . Long dashed lines: Flow with a more recent initial time 
ti = 10 Gyrs; Solid lines: Flow with stellar M/Lb increased to 16; Short dashed lines: Flow with a* reduced by 2 with M/Lb = 8. 

Fig. 4. — A comparison of single temperature (IT) spectra (for variable q M87 model) with a composite spectra of a single temperature 
plus cooling flow (IT -|- OF) with equal normalizations (for q = 1 M87 model). The OF spectrum extends from down to 0.1 keV. Upper 
panel: Spectra at il 30 kpc: (IT + CF) composite spectrum for gas at T/j = 3 keV (solid line); IT spectrum for gas at average temperature 
T = 2.5 keV (dotted line). Lower panel: il ss 1 kpc: (IT + CF) composite spectrum for gas at = 0.95 keV (solid line); IT spectrum for 
gas at T = 0.8 keV (dotted line). 

Fig. 5. — Gas density and temperature profiles for M87/Virgo at time tn including thermal conductivity reduced at all radii from the 
Spitzer or saturated value by a factor /. Upper panel: The effective density profile n(r) corresponding to emission from background and 
cooling gas. Variable q flow for / = 0.02: {heavy solid line) and / = 0.05 {heavy dashed line); Flow with q = 1 for / = 0.08: {light solid line) 
and / = 0.15 {light dashed line); Lower panel: Projected temperature T{R) profiles; line types are identical to those in the upper panel. 

Fig. 6. — Gas density and temperature in the q = I model for NGC4874/Coma at time t„. Upper panel: Density profiles n(r) with solid 
curves labeled as in Fig. 1. Long dashed line: /3-model gas profile from Briel et al. (1992) appropriate for large r. Lower panel: Temperature 
profiles for Th{r), T{r) and T{R) labeled as in Fig. 1. Data are from Arnaud et al. (2001) {open circles) and Vikhlinin et al. (2001) {filled 
circles). 

Fig. 7. — X-ray surface brightness distribution T,x{R) (0.5 - 2 keV) at tn in NGC4784/Coma before {solid lines) and after convolution 
with Chandra PSF {dashed lines). Upper panel: For q = 1 flow; Lower panel: For variable q flow. Circles are Chandra data from Vikhlinin et 
al. (2001) and squares are ROSAT data from Dow & White (1995). 

Fig. 8. — Gas density and temperature in the variable q model for NGC4874/Coma at time t„. Line types and data are identical to those 
in Figure 6. 

Fk;. 9. — Mass flows in M87. Variation of M{r) = fyaiTrr"^ in computed models {solid lines), Mi2{r) from Equation (12) {dashed lines) and 
integrated rate of cooling dropout M,-[a{r) {dotted lines), all evaluated at time tn = 13 Gyrs. Upper panel: For the g = 1 solution. Center 
panel: For the variable q solution. M (r) computed from the model (solid line) changes sign and becomes positive (outflow) for r > 53 kpc. 
Lower panel: An alternative variable q solution at tn in which M87 is evolved separately, then inserted into Virgo at t, = 6 Gyrs. 

Fig. 10. — Radial variation of M{r) {solid lines), Mj^a{r) {dotted lines) and Mi2{r) {dashed lines) for flows with thermal conductivity 
reduced by a factor /. Upper panel: q = 1 flow with / = 0.02; Lower panel: variable q flow with / = 0.08. 
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